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Abstract: One way to model the strange metal phase of certain materials is via a holographic 
^^ \ description in terms of probe D-branes in a Lifshitz spacetime, characterised by a dynamical 

f^ ' exponent z. The background geometry is dual to a strongly-interacting quantum critical 

^^ ■ theory while the probe D-branes are dual to a finite density of charge carriers that can exhibit 

t*^ ' the characteristic properties of strange metals. We compute holographically the low-frequency 

[~^ . and low-momentum form of the charge density and current retarded Green's functions in 

^^ . these systems for massless charge carriers. The results reveal a quasi-particle excitation when 

z < 2, which in analogy with Landau Fermi liquids we call zero sound. The real part of 
the dispersion relation depends on momentum k linearly, while the imaginary part goes as 
k> I fc^'^. When z > 2 the zero sound is not a well-defined quasi-particle. We also compute 

;h ' the frequency-dependent conductivity in arbitrary spacetime dimensions. Using that as a 



measure of the charge current spectral function, we find that the zero sound appears only 
when the spectral function consists of a single delta function at zero frequency. 

Keywords: AdS/CFT correspondence. Gauge/gravity correspondence. 



* E-mail address 
'E-mail address 
'■E-mail address 



choyosSphys . Washington . edu 
cQiobSmppmu . mpg . de 
jbnwuOitp . unibe . ch 



Contents 

|l|. Introduction and Summary |^ 

^. Probe Thermodynamics ^ 

^. Green's Functions at Low Frequency ^ 



Holographic Zero Sound |12 



gl| 1 < 2 < 2 |T2 



4.2 z>2 13 



|. AC Conductivity |4 

^. Discussion and Conclusion 15 



1. Introduction and Summary 

Gauge-gravity duality [1-3] provides a novel way to compute observables in strongly-coupled, 
scale-invariant systems at finite density, and thus may be useful for studying condensed matter 
systems near quantum critical points, such as the "strange metal" phase of some heavy fermion 
compounds and possibly also of some high- Tc materials. In particular, gauge-gravity duality 
may be able to reproduce the simple properties of strange metals, such as electrical resistivity 
scaling linearly with temperature T, which cannot be derived in a Fermi liquid description, 
which assumes the appropriate degrees of freedom are weakly- interacting quasi-particles. Of 
course the ultimate goal is to understand why some strange metals have a high Tc- 

Gauge-gravity duality is holographic: it equates a weakly-coupled theory of gravity on 
some spacetime with a strongly-coupled field theory living on the boundary of that spacetime. 
The spacetime symmetries of the field theory are dual to the isometrics of the bulk metric. 
Here we will focus on field theories invariant under the Lifshitz group, which includes scale 
transformations of the time coordinate t and spatial coordinates x of the form 

t^XH, x^Ax, (1.1) 

where A is real and positive and z is called the "dynamical exponent." The dual gravity 
theory then lives in Lifshitz spacetime, with metric [4] 

as = Qav dx^dx = -^ ^ H 5- , (1.2) 



where r is the holographic radial coordinate: r — t- oo corresponds to the infrared (IR) of 
the field theory, while the near-boundary region r — t- corresponds to the ultraviolet (UV). 
Notice that when z = 1, the metric becomes that of anti-de Sitter space (AdS), and the dual 
field theory symmetry is enhanced to the relativistic conformal group. 

One proposal for holographic model building of strange metals is to use probe D-branes in 
a Lifshitz spacetime [5]. In this approach, an appropriate configuration of fields on the probe 
D-brane represents a finite density of massive charge carriers [6-8], while the background 
spacetime represents a strongly-coupled neutral quantum critical theory. Holographic calcu- 
lations then reveal that such systems can indeed reproduce the simple properties of strange 
metals, including electrical resistivity linear in temperature [5]. However, ingredients beyond 
Lifshitz spacetime, for example a nontrivial dilaton, are needed to reproduce all the measured 
strange metal properties [5,9]. 

When z = 1 and T = 0, holographic calculations also reveal the existence of a propagating 
quasi-particle excitation producing a pole in the charge density and current retarded two 
point functions, which we denote G^(a;,A;) and Gj^{uj,k), where uj is frequency and k is 
momentum [10-12]. In a weakly-coupled Fermi liquid, fluctuations in the shape of the Fermi 
surface produce a collective excitation at zero temperature, Landau's so-called "zero sound." 
The dispersion relation of the mode in these holographic systems was identical in form to that 
of Landau's zero sound, with real part linear in k and imaginary part going as A;^, hence the 
mode was dubbed zero sound by analogy. We will call this mode "holographic zero sound." 
The physical origin of this mode remains mysterious since it appears in systems with both 
fermions and scalars — systems that do not fit neatly into either Fermi liquid or quantum Bose 
liquid theory. Indeed, given the existence of the holographic zero sound, as well as unusual 
scaling of the heat capacity with temperature, these holographic systems may be new kinds 
of (strongly-coupled) quantum liquids [10]. 

In this paper we explore the fate of holographic zero sound in these systems when z > 1. 
To do so, we compute holographically - and completely analytically, without numerics - 
the low-frequency and low-momentum forms of G^^{uj,k) and G^{uj,k), for massless charge 
carriers. When z ^ 2 we find 



w^ ^M . , s cj"^ 



aw(.,k) = ^G%i.,k)^^,-^-,—^, (1,3) 

where v and c are dimensionful constants that depend on the density.^ We obtain the 

holographic zero sound dispersion relation by setting the denominator of G**(a;,A;) to zero. 

Whether the mode is a quasi-particle depends on whether the cj^ or w^/^+^ term is larger at 

low frequency. When z < 2, the w^ term dominates, and the mode remains a quasi-particle, 

although with a dispersion relation modified relative to the z = 1 case: the real part remains 

linear, but the imaginary part goes as k"^'^. When z > 2 the real and imaginary parts are of 

the same order, and hence the mode is no longer a sharply-defined quasi-particle. 

^When z = 2 the last term in the denominator becomes CLo'^log{auj'^), with a scale a depending on the 
density. 



We also compute the AC conductivity associated with the charge carriers, generahzing the 
results of ref. [5], for two spatial field theory directions, to an arbitrary number of spatial 
directions. When z < 2, the AC conductivity's dependence on a; is a pole formally identical 
to the high-frequency or "collisionless" limit of the standard Drude conductivity, iu^^ [5], 
indicating that the charge current spectral function, which is a measure of the density of 
charged states, consists only of a delta-function term at a; = 0. When z > 2, the w-dependence 
is instead a power law of the form w"^'^, indicating a nontrivial density of states at low 
frequency. 

In short, we find that holographic zero sound appears in the absence of low-frequency 
charged states. That alone is perhaps not surprising: a continuum of states can "smear out" 
a pole. Indeed, in interacting systems with a Fermi surface, the zero sound pole is most 
easily identified if the zero sound velocity is greater than the Fermi velocity, so that the pole 
is outside a continuum of particle-hole states.'^ The change in behavior of the dispersion 
relation at z = 2 is the novel feature of the holographic systems. 

This paper is organized as follows. We review the thermodynamics of these systems in 
section ^. We compute the Green's functions in section |3|, and use the results to study the 
holographic zero sound and AC conductivity in sections ^ and ^, respectively. We conclude 
with some discussion and suggestions for future research in section p. 

2. Probe Thermodynamics 

In this section we review the results of refs. [8,10,15-17] for the thermodynamics of the charge 
carriers in these holographic systems. 

In gauge- gravity duality, a theory of gravity is dual to some large- A^, strongly-coupled 
non-Abelian gauge theory. Nf probe D-branes are dual to Nj fields in the fundamental 
representation of the gauge group(s) in the probe limit Nf ^ N [18]. We will call these fields 
quarks or flavor fields, in analogy with Quantum Chromodynamics. 

We will consider Nf coincident D-branes, whose action is the non-Abelian Dirac-Born- 
Infeld (DBI) action, describing the dynamics of the worldvolume scalars and U{Nf) gauge 
fields.^ The field theory then has a global U{Nf) flavor symmetry. We produce a finite 
density of charge carriers by introducing a chemical potential for the diagonal U{1) C U{Nf). 
The density operator J* is the time component of the conserved f/(l) current J^, and is 
holographically dual to At, the time component of the U{1) gauge field living on the D- 
branes. We will thus introduce in the bulk D-branes with only At{r) present. The action 
then reduces to the Abelian DBI action. 



^For textbook reviews of zero sound in Fermi liquids, see for example section 1.7 of ref. [13] or section 5.4 
of ref. [14] 

'^In string theory D-brane actions also include Wess-Zumino terms, describing the coupling of D-brane fields 
to background Ramond-Ramond fields. To keep the discussion general without committing to a specific string 
theory realization of Lifshitz spacetime [5], we will omit Wess-Zumino terms. 



We will only consider massless flavor fields. A flavor field mass operator would be dual to 
a scalar field on the D-branes. We omit any such scalar from our calculations, although we 
comment on massive flavors in section |^. 

We assume our D-branes are extended along q spatial dimensions of the Lifshitz spacetime.^ 
The Abelian DBI action for our D-branes is then 

S = -NfTovfclrdtd'^x ^-det [g^f, + {2TTa')Fab] , (2.1) 

where To is the tension, gab the induced metric, and Fab the U{1) field strength of the D- 
branes. The factor V is the volume of any internal space that the D-brane may be wrapping, 
and 2Tra' is the inverse string tension, which will typically be present in an actual string 
theory system. Inserting our ansatz At{r), performing the trivial integrations over dt and 
d'^x, and dividing out the subsequent (infinite) volume factors, we obtain the action density 
S (which for brevity we will call the action henceforth). 



S 



-MJdr 5^/2 ^|5ti|9rr-i;' , N = NjTdV , (2.2) 



where a tilde denotes a factor of 27ra', so At = (27ra')^t, and the prime denotes dr- The 
charge density in the field theory is 



(J') = %=M,lL-^M}£=. (2.3) 

Solving for A[{r), we find 



* \/|5it|5rr-A 



where d = {J^) /M . Inserting the solution for Aj(r) back into the action, we find 

S = -M[dr\gu\'/Vd'^4^- (2.5) 

J ygix + d^ 

We obtain the field theory chemical potential and free energy by performing the integrals 
for A!^{r) and S. The chemical potential is 

^'~Jo "^^ ^'^"^ - 2na' Jo ^' V^^^^T^ " 4vr«'g T (i) ' ^^'^^ 

The grand canonical (Gibbs) free energy density is 

/•oo -{2q+z+l) _ J 

n^ = _S = M dr =-M ^^. (2.7) 

Jo Vr-2g -\-(f z + q 



*If q is less than the total number of spatial dimensions, then in the field theory the flavor fields propagate 
only along some g-dimensional defect. 



These results are deceptively simple. We have actually removed various power-law divergences 
at the r = endpoint of the integrations via analytic continuation. To see the divergences 
explicitly, we expand the integrand of eq. (^J) ^°^ small r, 



/ dr ^ = dr r-(«+"+i) f 1 - -dV^" + O {d^r^") J . (2.8) 

The first term gives a density-independent divergence for any z, coming simply from the 
volume of Lifshitz spacetime. As we increase z to q and beyond, new density-dependent UV 
divergences appear. These are proportional to even powers of d, and appear to be associated 
with operators ( J*)^" with n G N that become relevant for larger values of z.^ We could have 
cancelled all of these divergences by introducing a regulator by integrating only to a cutoff 
r = e and then adding counterterms^ at r = e [19-21], but these counterterms generally 
also introduce additional finite density-dependent terms. Fixing these finite terms actually 
determines the thermodynamic potential, and thus the thermodynamic ensemble. We have 
implicitly fixed these terms via analytic continuation, and hence we have chosen an ensemble. 
A straightforward calculation reveals that for us 

(j') = -f^, (2.9) 

duo 

which shows that we are working in the grand canonical ensemble. Notice that with massless 
flavor fields, any finite chemical potential will produce a finite density. 

Analytic continuation cannot remove logarithmic divergences, which is indeed the case in 



eqs. ( pl^ ) and (2^) when z = (2n — l)q. These divergences occur because the operators 
(J*)^" are classically marginal at these z values [5], and if added to the action, can produce 
a breaking of scale invariance through quantum corrections, similarly to what happens for 
multi-trace scalar operators [22]. 
The energy density, Eq, is 

eo = no + fio{J') = -^no. (2.10) 

In the grand canonical ensemble and at zero temperature, the pressure Pq is simply Pq = 
— Oo = iz/q)eo, hence (dP/de)^ = z/q [12,17,23], a result fixed completely by Lifshitz scale 
invariance.^ However, this quantity is only the speed of normal/first sound in the relativistic 
z = 1 case. In general, the speed of normal/first sound is the derivative of the pressure 
with respect to the mass density, and the latter only coincides with the energy density in a 
relativistic system. In a system with z > 1, a velocity v must be dimensionful: given [oj] = z 



''The scaling dimension of the current operator is [J*] = q, while the volume element has dimension 
[dtd''x] — —z — q, so as we increase z to q and beyond, the operator J*,P goes from irrelevant to marginal (at 
z = q) to relevant. Similarly, when z = 3g, J*J^J*J* becomes marginal, and so on. Only even powers of J* 
can appear due to charge conjugation invariance J* — >■ —J*. 

^We write the counterterms explicitly in section H. 

'^Holographic calculations of (9P/9e)^ using asymptotically Lifshitz black holes yield the same result (see 
for example ref. [24]). 



and [k] = 1, we have [v] = [co/k] = z — 1. Indeed, we will later find a holographic zero sound 
speed proportional to S^^^'^''^. 

To study thermodynamics at finite temperature, we need a black hole in the bulk whose 
Hawking temperature is dual to the field theory temperature [25]. A number of asymptotically 
Lifshitz black holes have been studied recently - for a sampling, see refs. [24,26-35]. To 
capture the basic qualitative physics without committing to a specific system, we will again 
follow ref. [5] and write a generic Lifshitz black hole metric 

, 9 dr^ f(r) , 9 dx"^ 

d^" = Jl^2 - ^dt' + 7^ • (2-11) 

We do not need an exact form of /(r). We only need to know that /(r) has a simple zero at 
the horizon r = rn, and that the temperature is fixed by the regularity of the Wick-rotated 
metric at the horizon, 

/(r) ^ cofl-— V T = co^, (2.12) 

r-^rn \ Th J 47r 

where cq is a dimensionless constant that depends on the specific system. 
The chemical potential at finite temperature is 

j-TH J.-Z (\ z z r^^^\ 



-2q / -iq' 

2Tia'z \ 2{2q + z) d? \ d^ 



^o-:i^ 1-^,7:^:^^+0 r^ 11, (2.13) 



where the second line is a small-temperature expansion. The probe flavor contribution to the 
finite-temperature free energy is 




. r zd . . AirdT ^ I ^ 1 , ^ 

-AA^/i-AA , ^ , +0{ —^ , 2.14 



where in the second line we have used eq. (|2.7|), in the third line eq. (p. 13), and in the 



fourth line eq. ( ^.12]) . The entropy density at zero temperature is proportional to the charge 
density [10,15,17], 

^o = - h^ =^-—-^d. (2.15) 



Indeed, in a string theory system, this entropy should be equivalent to the entropy of a 
single quark, represented in the bulk by a single string, times the density d of quarks [15]. 
A nonzero entropy at zero temperature, or equivalently a degeneracy of states, indicates a 
potential instability [36], since typically any perturbation will break the degeneracy and the 
system will flow to a new, generically non-degenerate state. In other words, we may be 
working with states that are not the genuine ground state. The z = 1 cases appear to be 
locally thermodynamically stable, however [16]. The specific heat at finite temperature scales 
with T as 



cv 



^ ds 
dT 



2q 



d 



A linear scaling, as in a Fermi liquid, occurs when z 
Bose liquid, occurs when z = 2 [17]. 



(2.16) 
2q, while a scaling T'', as in a quantum 



3. Green's Functions at Low Frequency 

Zero sound should appear at zero temperature as a pole in the Fourier-transformed retarded 
two-point function of the density operator, G^ (to, k) [10]. In what follows we will also compute 
the density-current and current-current retarded two-point functions G*^{u}, k) and G^{uj, k). 
The latter will give us the AC conductivity. 

In gauge-gravity duality the on-shell action is the field theory generating functional [2,3], 
hence we need to take two functional derivatives of the on-shell bulk action, which in turn 
means we need to solve the linearized equations of motion for the fields dual to the density 
and current operators, namely the U{1) gauge fields At and A^ on the probe D-branes. The 
holographic zero sound pole will appear as a quasi-normal frequency of these fields. We thus 
consider fiuctuations a^ of the gauge fields, 



Mr) 



A^{r) + a^{r,t,x) 



(3.1) 



where because the field theory is invariant under rotations, we only need to consider fluctua- 
tions at and ax with r, t, and x dependence. We work in a gauge where ar = 0. To quadratic 
order in the fluctuations, the action is 



Sa^ 



M 



dr dt dx ql!} 



■?2 I |~'2 
9rrhx ~ \9tt\0.x 



+ 



gtt\grr — Af- 



\gtt\grrat 



gtt\gr 



K' 



3/2 



(3.2) 



where ftx = ^tOx ~ QxO't- Performing a Fourier transform. 



af,{r,t,x) 



duj dq 

(2^ 



-iwt+ikx 



a^{r,uj,k), 



(3.3) 



the linearized equations of motion become 
dr 



9/21 I / 

9 XX \gtt\grrO-t 



dr 



mWrr - A ) 
9/2-1 1 I / 

dxx \9tt\ax 



I a'2 

9tt\grr — ^t 



q/2-1 
Qxx Srr 



9tt\gr 



A 



+ 



q/2-1 
Sxx 9rr 



[k'^at + ojkax) = , 



[uj ax + ojkat) = . 



(3.4) 



(3.5) 



\gtt\gTT — A 

Additionally, we find a constraint (aj.'s equation of motion, written in a,. = gauge) 

grrgxxi^a't + ( \gtt\grr - A[^] ka'x = 0. (3.6) 



If we solve the constraint equation for a'^ and plug into the Ux equation of motion, we find 
the at equation of motion. We thus need only solve the constraint equation and at equation. 
We will work with the gauge-invariant electric field (in Fourier space) 



E{r, oj, k) = ojax + kat ■ 
The constraint equation and definition of E give us 
, u{r)'^k 



u{r)'^k'^ — uj"^ 



E', 



u{r)'^k^ — oj"^ 



E'. 



where 



/ N \9tt\grr-A[^ \gtt\ /, , ,2 -q\-l/2 

V 9rrgxx V gxx 



^{-l)(l + d2^2,^-l/2 



Plugging into the at equation, we find (after some algebra) an equation for E, 



E" + 



9j.ln 



' (g-3)/2 -1/2, I 

gxx grr \gtt\ 
u {k'^v? — oj'^) 



m\ 



In terms of E the quadratic action is 

M r (<?-3)/2 1/2 

bg2 = — ar auj ak 

2 ./ u 



E'^ + 



\gtt\ 



-E 



Qrriu'^k'^ — W^) 

Introducing a cutoff at r = e and integrating by parts, we obtain in the e — )■ limit 



Sa'^ 



M 



cZ-q+l 



, duj dk — -— - 
2 / A;2 



-E{e)E'{e). 



(3.7) 



(3.8) 



(3.9) 



(3.10) 



(3.11) 



(3.12) 



We now need to solve the equation of motion with an ingoing boundary condition im- 
posed at the r — )■ 00 "horizon" [37-39], insert the solution into 5*^2, and then functionally 
differentiate to obtain the retarded correlators, for example 



Gl{u:,k) 



6at{e 



— *? 

:2 "■ 



5E{e)V 52 



6at{e)J 6E{e 



2^a 



(3.13) 



Defining 

i 

6E{e 
we get 



n(^, k) = T^^n^S,^ ' (3-14) 



G*}^{uj,k) = k^Il{uj,k), G''^{oj,k)=ujkU{uj,k), G^j^ {oj , k) = uj^ U{u} , k) . (3.15) 

Our goal is to obtain (the low-frequency and low-momentum form of) the variational deriva- 
tive n(a;, k), which determines the retarded two-point functions that we want. 

We have not been able to solve eq. ( |3.10| ) exactly for all values of uj and k. A numerical 
solution is of course possible, however we can obtain the low-frequency behavior of Hiuj, k) 
by solving eq. (|3.10| ) in two different limits and then matching the two solutions in a regime 
where the limits overlap, following ref. [10]. Specifically, we first obtain a solution by taking 
r large and then taking the small frequency and momentum limit, which means wr^ <^ 1 and 
kr <^ 1 with ujk~^ fixed. We then repeat the process, now taking the small frequency and 
momentum limit first followed by taking r large. 

Taking r — t- oo, the equation for E becomes 

E" + -\2-(z- 1)1 E' + ^^2(^-1)^ = . (3.16) 

r 

The solution we want involves a Hankel function of the first kind, 

which for large r describes an in-going wave: E ~ r~^ ^luir^/z g^g ^ ^, ^^^ Now we want to 
take frequency and momentum small, which means in particular cor^ <C 1, since the solution 
above does not depend on k. The form of E in this limit depends on the value of z. When 
z ^ 2, E approaches 

..cr(l.l)-'(,-..a„9-£r(l-l)(|i)'-.-. ,3.., 

When z = 2, the wr^ <^ 1 limit of E involves a logarithm 

E^C + C—\log(ujr^)-log4 + -f], (3.19) 

where 7 is the Euler-Mascheroni number. 

Next, we start by taking uir^ <C 1 and kr <^ 1 with ujk~^ fixed. The equation for E 

becomes 

'2-q-2{z-l) u' Sk'^u^ - uj"^' 

r u k'^v? — bj^ 



E" + 



E' = 0. (3.20) 



The solution is 



E = Ci + C2{k' 



^q-z 



-2+q+z 



-UJ 



-2 + q + z 



■ 2F1 



2FI 



Q-z 3 q 

[ 2q '2'' + 

"-2 + g + z 1 


— z 
2q 

1 + 


2q '2 



;-dV'' 



-2 + q + z^ 
2q '■' 



2^2g 



-d^r 



(3.21) 



where Ci and C2 are constants. We can match to eqs. (3.18) and ( 3.19| ) by taking r large. 
The form of E in this limit also depends on the value of z. When z ^ 2, the large r limit of 
E is {B{x, y) is the Beta function) 



E^Ci + 



9i 

d 



zk di 



B 



z z 
2^'2^ 



w^d-J^/l z-2 2-z 

B\ - H , 

2q \2 2q 2q 



+ 



uP'r^ 2 



(2-z) 



(3.22) 



which, as a function of r, consists of a constant and an r^ ■^ term, which indeed match the 
terms in eq. ( 3.1^ ). When z = 2, the large r limit of E changes only in the w^ term. 



E«Ci + 



9l 
d 



ed^/i 



1 1 1 



UJ 



B{---,-] log(2dr 



(3.23) 



2 q qj q 

which again consists of a constant and a logarithm, as in eq. ( |3.19| ). To do the matching, we 
rewrite the logarithmic term as 

2d 



log(2dr^) = log(^(a;r2^«/' 



UJ 



q/2 



log 



,./2;+|l°s^--'^ 



(3.24) 



The boundary behavior of the solutions is given by the small r expansion of the solution in 
eq. ( [3.211 ): E ~ Ci + k'^C2r''~^ /{q — z). When z < q,we have that to leading order E{€) ~ Ci, 
so 



n(a;, k) = lim ■ 



. >'6 <5^(e)2 ^"^ 
while to leading order E'{e) ~ k'^€'^^^^^C2, so*^ 



Sa2 



TV f e^~'^+^ - 

-— dojdk—^E{e)E'{e) 






e->0 



{2TTaY dujdkCiC2. 



(3.25) 



(3.28) 



*When z > q, C2 becomes the source in the generating functional, and, as mentioned in section Q, 8^2 has 
new divergences as e — >■ 0. To make 8^2 finite we add a boundary term, which for z > q is 



_ fdud.^^ 



tt XX f^2 ( \ 

'111 E (e) 



N 



dwdk 



Q- ^^z-9g,2 



fc2 



and for z — q, where E{e) ~ Ci + C2fc^ log e, is 



S, 



M 



dujdK- 



1 



-l''l-E\.)^^jd^dn^^^^^^ 



E\e), 



' E\e), 



(3.26) 



(3.27) 



2 J " ' k^loge 
where 7^;^ is the induced metric a,t r — e. Adding such a term is equivalent to a Legendre transformation. 



which makes Ci the source again. The result for 5^2 is then identical to eq. (3.25), except with opposite sign. 
The boundary term does not change i?'s equation of motion, so the solutions and matching are the same, while 
the overall sign is not crucial for our physical arguments. The z > q cases are thus identical (up to a sign) to 
the z < q cases, so we can freely take 2 > q in our results. 



10 



To find Il{u},k), we need to find C2 in terms of Ci, which can be done by matching the 
constant term and the coefficient of the r^~'^ term in eqs. ( |3.18 ) and ( p.22 ). More specifically, 



the matching gives us two equations involving C, Ci, and C2, and we use one equation to 
eliminate C so that the remaining equation gives us C2 in terms of Ci. When z = 2, we 
match the constant term and the coefficient of the logarithmic term log(a;r^) in eqs. ( p.l9|) 
and (3.23|). After that we immediately obtain Il{uj,k). 



Assembling all the ingredients, we have our main result, the low- frequency and low- 
momentum form of n(a;, k): 

U{u;,k)^—-^ ^ pTTT' (3-29) 

where 

GoM= ,, , ., ^„ , (3.30) 



w log (aw^) 



and oi, 02 and 03 are constants. For z ^ 2 they are 



2-z /I 1\ /I 1\ dV /I z-2 2 

«2= . .. 0/ r --- r - + - -— ij - + 



7r(2z)^^^/'' \z 2y V-2 2y 2g V2 2g ' 2g 

vr 
as = i + tan — , (3.31) 

where for later reference we note that ai and a2 are real while as is complex. For z = 2 the 

constants are 

4^2/9 /I 1 1\ 1 
"1 = -B ,- , a2 = t, 03 = , (3.32) 



^ q qj vr 

and another (dimensionful) constant appears inside the logarithm in Gq{uj), 

a = ^{2d)-^/'^ . (3.33) 

In the denominator of Ii{u},k), when z < 2 the w^ term is larger than the Gq{uj) term, 
whereas when z > 2 the Gq{u}) term is larger. The dominant term dictates the low- frequency 
behavior of the retarded Green's functions, and hence of the holographic zero sound and AC 
conductivity, as we will show in detail in what follows. 

Notice that z = 2 is a special value regardless of the number of spatial dimensions q, which 
is similar to what is observed for the worldvolume scalar that describes the embedding of the 
D-brane (without a magnetic field) [40]. To see more clearly how this happens in the bulk. 



let us consider the near boundary limit r — )• of eq. ( 3.20 ): 



E" + i±^^^' = . (3.34) 
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We can compare this to the equation near the horizon, eq. ( 3.16 ), when a; — )• 0: 



E" + ^-^E' = 0. (3.35) 

r 

The near-horizon equation looks hke the near-boundary equation with an effective number of 

spatial dimensions q ^ q = 2(z — 1). This suggests that the finite density theory flows towards 

an IR fixed point that depends only on the value of z, much like the holographic system of 

ref. [41], but restricted to the probe D-brane sector. If we also identify q with the effective 

IR dimension of the density operator, then the operator ( J*)^ becomes marginal when q = z, 

which happens precisely when z = 2. More explicitly, when z < 2 and q < z, ( J*)^ is relevant, 

while when z > 2 and q > z, (J*)^ is irrelevant, so z = 2 would separate two different kinds 

of IR effective theories. Similarly to the logarithmic UV divergences we encountered in sec. y, 

the appearance of the logarithm in Go{io), which breaks scale invariance, may be due to (J*)^ 

becoming marginal at z = 2 in the IR effective theory. 

4. Holographic Zero Sound 

The denominator of 11(6^;, A;) vanishes - and hence the retarded Green's functions have a pole 

- whenever 

k{u)) = ±^=^a2UJ^ + asGoiuj) . (4.1) 

vol 



This determines the dispersion relation of the holographic zero sound mode: we need only 
expand the right-hand side of eq. ( [4.1[) and then invert to find u!{k). As mentioned above, 
the value of z determines whether the co'^ term or Gq{uj) term dominates at low frequencies, 
so we will consider two cases in turn. 

4.1 1 < z < 2 

For 1 < z < 2, the uj'^ term under the square root has larger norm, so we expand k{uj) as 



V ai V «2 



= ±iO 

and then invert to find 



1 + ^^-1+2^ + 0^-2+4/. 

2a2 V 



l/z 



u{k) = ±k,[^-^(^\ ' k^/^ + O (k-'^'/A . (4.2) 

The mode behaves as a quasiparticle since the imaginary part goes as k"^'^ , which is smaller 
than the real part at low momentum. As a check, if we set z = 1 we recover the result of 
ref. [10], 

u.{k) = ± A _ , ,f"^^\. k' + O {k') . (4.3) 

^ ra-4 r^ 
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The speed of the holographic zero sound, vq, is given by 

As promised in section |2|, when z > 1 the speed of the holographic zero sound is dimensionful. 
In the relativistic case, z = 1, the speed of holographic zero sound coincides with the speed 
of normal/first sound. For 1 < z < 2, the poles iii Tl i — ^ J and Tl ^-^ J determine whether 

Vq remains finite or goes to zero as z — ?• 2, with the details depending on the dimension q. 
When q = 2, the F functions (and hence their poles) immediately cancel and Vq = |d^~^, 

which is clearly finite as 2: — )• 2. For any 9 > 2, F( 2 — ^ ) remains finite but Fl -^ 1 has a 
pole as 2; — >• 2, so Vq goes to zero from above. 

4.2 z> 2 

We consider first z > 2. Now the Gq^oj) term under the square root in eq. ( [4.1D has larger 
norm, so we expand k{uj) as 

k{uj) = ± — w 2z 1 -I UJ z 

V"l/ V "3 / 

(4.5) 



1 
a3 \ 2 z+2 



±1 — 1 LO^ 
"1 



"3 



and invert to find 



z 2(z-l) 

Note that the leading term has a coefficient that is complex. In particular, the real and 
imaginary parts are of the same order, hence the excitation is not a quasi-particle. 
Finally, when z = 2 we have 

k = ±^^ v^a2 + as log ( a uj'^). (4.7) 

vol 

Expanding for small uj we find 

kicj) = ±^ V«3log(aw2) _ _^^ (03 log (aa;2))"^ +oL log"! (ai^^^) . (4.8) 



Although this mode is gapless, the dispersion relation differs from the holographic zero sound 
mode by logarithmic factors. 
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5. AC Conductivity 



In this section we provide some physical intuition for the change in behavior of the holographic 
zero sound mode as a function of z by studying the low-frequency AC conductivity (t((jj), 
defined as 

(T(a;) = --G^^(L^,/fc = 0). (5.1) 

Recalling that the imaginary part of the retarded Green's function is proportional to the 
spectral function of J^, we see that the real part of the conductivity provides a measure of 
the density of states that couple to J^ . 

To obtain aiui) for our system, we first set /c = in our result for n(w, k) in eq. ( |3.29| ) and 
then plug into eq. ( |3.15 ), which gives 



G^^(w,A: = 0)oc 



Mu? 



w^O 



a2ijJ^ + a3Go(a;) 



-M 



a. 



z <2 

aj^ (log(aa;2))"^ z = 2 



an UJ 



l.,,l-2A 



z>2 



(5.2) 



and then obtain the AC conductivity at small frequency from eq. (5.1) 



UJ 



ioo UJ 



z <2 



ia^, '^ (a;log(au;^)) z = 2 
z>2 



(5.3) 



ia^^oj-^/' 



Our result for cr{uj) not only generalizes that of ref. [5] to any number of spatial dimensions 
q, but is an independent confirmation of the result, since we derived it via a slightly different 
route from that of ref. [5] . 

Recall that when z ^ 2, 02 is a real number while 03 is complex. When z < 2, the 
conductivity is purely imaginary and has a simple pole at zero frequency. A Kramers-Kronig 
relation then implies that the real part of the conductivity, and hence the spectral function, 
consists only of a delta function at zero frequency. The system has no charged states at small 
frequencies. Notice that the iuj~^ behavior is expected in a system with free charge carriers 
(or with translation invariance), and indeed is formally identical to the high-frequency or 
"coUisionless" limit of the standard Drude result. When z > 2 the conductivity, and hence 
the spectral function, has a power-law dependence, and no clean holographic zero sound 
quasi-particle exists. 

Something analogous happens in interacting systems with a Fermi surface [13,14]. Here 
the charge density spectral function includes a continuum of particle-hole pairs at frequencies 
Lo < vpk, where vp is the Fermi velocity. The zero sound gives a contribution at uj = v^k. 
If vq > vp, then the zero sound is outside the continuum, and will produce a sharp quasi- 
particle peak in the spectral function. However, if vq < vp, the zero sound mode suffers 
Landau damping through scattering with particle-hole pairs, and hence will be reduced to 
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a broad resonance with no good quasi-particle interpretation. A similar mechanism appears 
to be at work in the holographic models we have studied: whenever the spectral function 
is non-vanishing at small frequencies, the holographic zero sound mode becomes a broad 
resonance. 

Notice that previous calculations of the DC conductivity from holographic probe brane 
systems yielded a finite result [7] , whereas the naive a; — )■ limit of our result for the AC con- 
ductivity is infinite for any value of z. This is not a contradiction. The different results come 
from the fact that the a; — )■ and T — )■ limits do not in general commute, or equivalently 
the ixjjT — 7- and bjjT — )■ oo limits are physically distinct [42]. The calculation of ref. [7] 
involved taking a; — )■ first, which is the so-called collision-dominated or hydrodynamic limit 
ujjT — 7- 0, while our calculation involves taking T — )■ first, which is the collisionless limit 
w/T —7- oo. The latter is the appropriate limit in which to see zero sound in Fermi liquid 
systems [13, 14]. 

6. Discussion and Conclusion 

We determined the low-frequency and low-momentum form of the density and current re- 
tarded Green's functions for massless charge carriers in the holographic model of strange 
metals proposed in ref. [5]. Our main physical results are that for z <2 the system's spec- 
trum includes a quasi-particle - the holographic zero sound - whose dispersion relation we 
computed, whereas when z > 2 the quasi-particle is "washed out" by a continuum of low- 
frequency states. In addition, we recovered the results of ref. [5] for the AC conductivity 
from this single contribution to the spectral function. We also speculated that the change of 
behavior at z = 2 may be due to the existence of an infrared fixed point where the effective 
dimension of the charge density is 2(z — 1), so that the operator (J*)^ becomes marginal at 
z = 2. 

Our z < 2 results are similar to those discovered in another holographic system, (3+1)- 
dimensional gravity in an extremal AdS-Reissner-Nordstrom geometry [43,44]. The dual field 
theory is a strongly-coupled (2+l)-dimensional CFT at finite density and zero temperature. 
Notice that here the density is not in the probe sector. Holographic calculations reveal a 
sound excitation in the spectrum, with a dispersion relation that is almost the T — )■ limit 
of the normal/first sound dispersion: the speed squared is still equal to the conformal value 
1/2, while the sound attention constant differs by about 10% [44]. The leading contribution 
to the low-frequency AC conductivity is a pole in the imaginary part [43]. 

The field theory states dual to extremal AdS-Reissner-Nordstrom and to our solutions share 
two important features. First, both have a nonzero entropy at zero temperature. Second, 
both appear to have an emergent scaling symmetry in the infrared, indicating some fixed 
point involving strongly-coupled degrees of freedom [41]. Notice however that holographic 
zero sound has been found in systems with zero entropy at zero temperature [23], which 
means that a nonzero T = entropy is not a necessary condition for the existence of a 
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holographic zero sound mode. The absence of a gap for charged excitations appears to be a 
necessary condition, however [11]. 

An alternative holographic approach to modeling strange metals is to introduce a charged 
probe fermion into AdS-Reissner-Nordstrom [41,45-47]. The probe fermion is dual to a 
charged fermionic operator whose two-point functions reveal a Fermi surface with either 
Fermi-liquid or non-Fermi-liquid properties, depending on the bulk fermion's mass and charge 
[41,45,46]. Given that the charge density is not in the probe sector, the fermions can only 
influence the charge and current Green's functions via back-reaction or one-loop effects. A 
calculation of the latter reveals that the probe fermions' contribution to the resistivity can 
scale linearly with temperature [47]. An open question is whether back-reaction or one-loop 
effects produce a zero sound pole, which may be a genuine Landau zero sound, corresponding 
to fluctuations in the shape of the Fermi surface. 

A related approach is to introduce a weak coupling between a Fermi liquid and a critical 
sector described by the near-horizon region of the extremal AdS-Reissner-Nordstrom black 
hole [48]. In this case the strange metal behavior comes from the critical sector's contribution 
to the fermion self-energy. A zero sound mode would receive corrections to its self-energy 
as well, modifying its dispersion relation. We have found indications that the near-horizon 
region of the probe D-branes also describes an IR critical theory, so we expect that an analysis 
similar to that of ref. [48] could be possible here, and that for the zero sound it would produce 
results similar to ours. 

Returning to the model of ref. [5], many avenues remain open for future research. As 
emphasized in ref. [5], a realistic analysis requires massive charge carriers, with a mass gap 
much larger than the temperature. In the bulk that means introducing a scalar field dual 
to the flavor mass operator on the probe D-brane worldvolume. In relativistic string theory 
systems, a special cancellation of metric factors occurs in the DBI action, allowing for an 
analytic solution not only for A[(r) but also for the scalar [8]. These analytic solutions 
then allow for an analytic calculation of the holographic zero sound dispersion relation [11]. 
The cancellation of warp factors does not occur for generic z, at least with a trivial dilaton. 
Nevertheless, we expect the physics for z > 1 to be similar to the z = 1 case, at least 
whenever a holographic zero sound quasi-particle exists. For example, we expect the speed 
of the holographic zero sound to go to zero as the chemical potential approaches the value of 
the mass from above [11]. 

In relativistic string theory settings, the holographic zero sound quasi-particle is absent at 
any finite temperature, apparently as a consequence of the probe limit [12]. In these systems, 
a finite temperature produces an energy density of order N'^, which is enormous compared to 
the probe flavor's order NfN energy density. The low-frequency dynamics is thus dominated 
by adjoint-sector physics, which induces charge diffusion. The spectral function grows linearly 
with frequency in this case, so the finite-temperature width of the holographic zero sound is 
big and the pole is smeared into a broad resonance. Presumably the same is true for any 
z < 2. A natural question is what happens if we leave the probe limit. Does the holographic 
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zero sound remain a sharp quasi-particle when Nf is on the order of A^? Notice also that a 
fully back-reacted calculation would unify several previous calculations, at least in systems 
with z = 1: at small field strength, where the DBI action is well-approximated by a Maxwell 
action, the results must reduce to those of ref. [44] (at least in (2+1) dimensions), while in a 
probe limit the results must reduce to those of ref. [10]. Notice that in both limits the speed 
of the sound mode was identical to the speed of normal/first sound. 

Lastly, the model building proposed in ref. [5] is worth pursuing (see for example ref. [9] ) . As 
noted there, some studies of the high- Tc cuprate strange metals find that the low-frequency 
AC conductivity scales as uj^'^''^, which suggests that any holographic model reproducing 
such scaling may not exhibit a zero sound pole. A holographic system that reproduces all 
measured strange metal properties will be more complicated than just probe D-branes in a 
Lifshitz geometry however, so the question of whether a holographic zero sound quasi-particle 
exists in these more complicated models should be studied on a case-by-case basis. 
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